Strong ExB shear flows in the transport barrier region in H mode plasma 
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We report the first experimental observation of stationary zonal flow in the transport barrier region 
of the H mode plasma. Strong peaks in Er shear mark the width of this region. Strong m = n = 
low-frequency (/ < 0.6 kHz) zonal flow is observed in regions of increased E^, suggesting substantial 
contribution of zonal flow to the spatial modulation of Er radial proflles. Radial localization of the 
zonal flow is correlated with a region of zero magnetic shear and low-order (7/5) rational surfaces. 
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Transport barriers (TBs) are radially localized regions 
in toroidal plasma where radial transport of particles or 
energy is drastically reduced. In the high confinement 
mode (H mode) [H, the presence of a TB is manifested 
as a steep density (or temperature) gradient near the 
plasma boundary. The top of this region is sometimes 
referred to as a pedestal. 

Characteristics of the H mode edge TBs are impor- 
tant. Spatial structure of a TB is closely related to the 
global stability, confinement and the plasma performance 
(for a review see, for example, 2]). Understanding and 
predicting characteristics of the TBs has become a focus 
of the international fusion community (see, e.g., 0, |3])- 
The ultimate goal of these studies is the optimization of 
the radial profiles of the plasma parameters in the future 
fusion reactor 5]. 

The formation of a TB has been ascribed to the gen- 
eration of a sheared radial electric field (or ExB flow, 
where E is the electric field and B is the magnetic field) 
which leads to the reduction in turbulence and transport 
0, 0, @1 ■ However, the physics of the TB formation is not 
yet well understood. Experimental studies of TBs are re- 
stricted due to difficulties in measuring radial parameter 
profiles with sufficient spatial and temporal resolution. 

In this Letter we report detailed experimental studies 
of the TB structure in H mode of the H-l heliac. It is 
shown for the first time, that distinct features in the elec- 
tron density profile, marking the pedestal and the foot 
of TB, spatially coincide with radially localized strongly 
sheared ExB flows. These radial regions are also iden- 
tified as regions where strong stationary m — n — zonal 
flows are localized in H mode. The radial localization 
of zonal flows also coincides with the position of a low- 
order rational surface and a minimum in the magnetic 
shear. These results confirm, to some extent, a hypoth- 
esis based on results of the gyrokinetic simulations that 



strong zonal flows developing near rational surfaces can 
provide a trigger for the TB formation 9]. 

We present results obtained in the H-l toroidal heliac 
[lo| (major radius of i? = 1 m and mean minor plasma 
radius of about (a) « 0.2 m) under the following plasma 
conditions (see, for example, and references therein): 
rie = 1 X 10^^ m-^, Te - 10 eV, Ti - 40 eV in argon at 
filling pressure of (1 — 4) x 10~^ Torr and at low magnetic 
fields, B — (0.05 — 0.12) T. Such plasma is produced by 
~ 80 kW of the radio-frequency waves at 7 MHz. Several 
combinations of Langmuir probes (single, triple probes) 
are used to characterize plasma parameters, such as the 
electron density, electron t emp erature, and electrostatic 
potential, as described in [ij]. Probes are also used to 
characterize poloidal and toroidal wave numbers of tur- 
bulent fluctuations UM- 
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The high confinement mode observed in H-l 1J| is sim- 
ilar to H mode in tokamaks. Typical electron density and 
plasma potential profiles are illustrated in Fig. [1] for low 
(L) and H modes. L and H modes are achieved respec- 
tively above and below critical magnetic field IJ] . When 
the magnetic field is close to the critical value, sponta- 
neous L-H transitions are observed. In these discharges, 
a triple probe is used to measure the electron density and 
potential on a shot-to-shot basis. Excellent reproducibil- 
ity of the measurements allows reliable determination of 
profiles without perturbing plasma by the probe arrays. 

Despite large differences in electron temperature, den- 
sity, and magnetic field, plasmas in H-l and in the TB re- 
gions of large tokamaks, are dimensionally similar. This 
dimensional similarity has been discussed in [isj . where 
it has been shown that the width of TB measured in ion 
gyroradii is very similar to that in, for example, DIII- 
D tokamak. However in absolute units, the TB width 
in H-l is substantially broader (30-40 mm) than that 
in larger experiments with stronger magnetic fields and 
lighter ions. This, in combination with low electron tem- 
perature in H-l, opens an opportunity to study structure 
of the TB using probes with sufficiently high spatial res- 
olution. 

The development of the TB in H mode plasma in H- 
1 is illustrated in Fig. [1] Radial profiles of the electron 
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FIG. 1: Radial profiles of (a) electron density, (b) plasma po- 
tential in L mode, and (c) electron density, (d) plasma poten- 
tial in H mode, respectively. Dashed guide lines and shading 
are used to mark two radial regions in the transport barrier 
(I and II). 
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FIG. 2: Radial profiles of (a) radial electric field and (d) shear 
in the radial electric field, computed using plasma potential 
profile of Fig. Wid). 



density and the plasma potential in L mode are rather 
featureless as seen in Fig. [TIa,b). In H mode, the cen- 
tral density doubles while the plasma potential becomes 
more negative in the central region and more positive at 
the edge (Fig. [Hc,d)). The increase in the density co- 
incides with the formation of the characteristic kink in 
the density profile at about p — r/a « 0.6, referred to as 
the pedestal. The rig profile outside the pedestal can be 
approximated by a straight line (Fig. (TJc)). The profile 
of the plasma potential (p also shows two characteristic 
kinks: one at the top of the TB, and the other at p = 0.8, 
which we will refer to as the foot of the TB. The third 
kink in the plasma potential is seen near the last closed 
flux surface (p = 1.0) and is due to the reversal of the 
radial electric field from negative (inside) to positive (out- 
side) . We use dashed guide lines and shading throughout 
the paper to mark two radial regions of interest: (I) - a 
region between the top and foot of the TB, and (II) - a 
region between the foot and the last closed flux surface. 

Profiles of the radial electric field and its shear E'^., 
derived from the radial profile of the plasma potential are 
shown in Fig. [21 Since is computed by differentiat- 
ing radial profile of the plasma potential, the (negative) 
maxima of the radial electric field can not be determined 
exactly. Three Er regions are seen: slightly positive Er 
inside the top of the transport barrier, substantial nega- 
tive Er ~ -1 kV/m in region I, and even more negative 
Er « -4 kV/m in region II. Correspondingly, the has 
distinct peak at the top, at the foot of the transport bar- 
rier and at the last closed flux surface. 

Fluctuations in the electron density and potential are 
strongly reduced in a broad range of frequencies from 
L to H mode as discussed in [16|. However, the low 



frequency (/ < 0.6 kHz) spectral feature increases in 
some radial regions, which will be discussed later. The 
power spectra of the fluctuations in the plasma potential, 
P{4>), at various radial positions in H mode are shown in 
Fig. El^a). The low-frequency feature (0.1 — 0.6 kHz) is 
dominant in H mode. Poloidal wave number kg of this low 
frequency component is measured using two poloidally 
separated probes. Measured poloidal wave number of 
ke = (2- 5) m-i at / = (0.1 - 0.6) kHz is indicative 
of the mode number to = 0. The toroidal mode num- 
ber is estimated using toroidally separated probes, as de- 
scribed in [3], and shows n — 0. Hence the strong low 
frequency fluctuations in the plasma potential are identi- 
fied as stationary zonal flows. It should be noted that it 
is usually difficult to align toroidally separated probes to 
exactly the same poloidal position. As a result, a phase 
shift between toroidally separated probes will occur due 
to the uncertainty in the poloidal separation between the 
probes. Ay: 



A^{f) = k«{f)ALn+ke{f)Ay, 



(1) 



where AL|| and Ay are toroidal and poloidal separa- 
tion between the probes respectively, and kg{f) is known 
from the phase difference between the probes which are 
poloidally separated. In case of a zonal flow, m = 0, the 
second term on the right-hand side becomes zero (since 
kg = 0), such that the poloidal uncertainty Ay becomes 
unimportant and the toroidal wave number can be reli- 
ably estimated by measuring Atp. 

Spectra similar to those in Fig. [Ha) have also been 
observed in the Compact Helical System (CHS) using 
heavy-ion-beam probe 17[. In that experiment, low fre- 
quency potential structures were also identified as sta- 
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FIG. 3: (a) Power spectra of the plasma potential in different 
radial regions; (b) radial profile of the spectral power density 
of stationary zonal flows (0.1 ^ 0.6 kHz). Hatched boxes 
indicate radial positions of the E,. maxima in regions I and 
II. 



tionary zonal flows. 

The spectral power density of the zonal flow (shaded 
spectral region of / = (0.1 — 0.6) kHz in Fig. [3fa)) varies 
along the radius. Radial proflle of the spectral power den- 
sity of the zonal flow in H mode is presented in Fig. [3Jb) . 
Stationary zonal flow is a band-like structure localized in 
radial region of 0.6 < p < 1.0. 

Two hatched boxes drawn in Fig. [3^b) indicate the 
uncertainty in the radial positions of the (negative) Ej. 
maxima in regions I and II ( Fig. [DJa)). It can be seen 
that the zonal flow maximum spatially coincides with the 
maximum in (negative) Er. This suggests that stationary 
zonal flow directly contributes to mean and may be 
responsible for the "corrugation" of the Er profile seen 
in Fig. Ha). 

The list of spatially coinciding phenomena in this 
plasma is complemented by the observation that the TB 
region appears in the vicinity of the zero magnetic shear 
in this magnetic configuration. The computed radial pro- 
file of the rotational transform t = 1/q (where q is the 
safety factor) is shown in Fig. [4l In addition to zero shear 
at p « 0.75, t — lA = n/m = 7/5 rational surfaces are 
present in both zones I (at p « 0.65) and II (at p ~ 
0.85). The accuracy of the t computation has been veri- 
fied using experimental electron beam mapping [l9| . The 
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FIG. 4: Radial profile of rotational transform in magnetic 
configuration discussed in this paper. 



existence of the 7/5 rational surfaces is also confirmed by 
the observation of the m — 5 chain of magnetic islands 
in the region of p « (0.83 — 0.87). The plasma current in 
the H-1 heliac is negligibly small (~ 10 A) and does not 
affect the vacuum magnetic structure. 

A possible role of low-order rational surfaces in the 
formation of H mode has been recognized since the first 
observation of H mode in stellarators [18, 20]. Spatial 
correlation of the rational surfaces with the stationary 
zonal flows, seen in Fig.[3ljb) and Fig.lH may be indicative 
of the generation of stationary zonal flows due to the 
influence of the rational surfaces, as suggested in 2l| . 

Formation of the TB and strong stationary zonal flow 
is also observed during spontaneous L-H transitions in 
H-1, described in 15|. Fig. [5l^a) shows temporal evolu- 



tion of the mean plasma density during spontaneous L-H 
transition. In this discharge the mean electron density 
jumps from about 0.6 x lO^'* m^"^ to almost 1.2 x 10^^ 
m~'^ in about one millisecond. Similarly to stationary H 
mode discharges described above, strong zonal flow in the 
TB region is observed in the H mode stage of discharges 
with spontaneous transitions. The spatial correlation of 
the TB regions and stationary zonal flow in the sponta- 
neous transitions is also observed. 

In Fig.[5Ub), the change of the fluctuation power spec- 
tra across the L-H transition at the radial position of 
p — 0.65 is illustrated. It can be seen that across the 
transition, fluctuations in the broad spectral range from 
0.6 to 100 kHz are reduced, while the spectral power 
of low-frequency zonal flow, / ~ (0.1 — 0.6) kHz, is in- 
creased. 

Many of the ingredients of the TB physics presented in 
this paper have been discussed with regard to H modes in 
tokamaks and stellarators. For example, the role of the 
low-order rational surfaces [1, , the role of zonal flows 
in L-H transition [2^, modification of the £V profiles by 
non-neoclassical (turbulence-driven) effects p, 24, 2^. 



Here we present for the first time experimental evi- 
dence that zonal flow, which develop near the n/m = 
7/5 rational surfaces, is spatially correlated with distinct 
regions of the radial electric field inside the TB region 
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wave is responsible for the sustainment of zonal flow in H 
mode, as has been proposed in [2^. Further theoretical 
work in this direction is needed. 

The authors would like to thank Santhosh Kumar for 
providing computed data on rotational transform. 



FIG. 5: (a) Time evolution of the mean plasma density dur- 
ing the spontaneous L-H transition, (b) power spectra of the 
plasma floating potential across the L (dashed line) to H (solid 
line) transition at p = 0.65. 



in the H mode plasma. The resulting strong peaks in 
Er shear coincide with the kinks in the density profile 
in H mode, suggesting that the strong peaks in the E'^. 
define the position and the width of the TB. Similar phys- 
ical picture has recently emerged as a result of analysis 
of gyrokinetic simulations of Dlll-D tokamak discharges 
Corrugations in the radial profiles of electron den- 
sity, temperature and radial electric field in DIII-D are 
observed near low-order rational surfaces. The develop- 
ment of strong zonal flows and strong Er shear layer in 
these plasma regions suggests the development of zonal 
flow as a trigger for the TB formation. 

Since zonal fiows are usually thought of as turbulence- 
driven flows, possible mechanisms of the zonal flow en- 
hancement and sustainment in H mode, when the level of 
turbulence is substantially reduced, need to be explained. 
It has been suggested in [l^, that the redistribution of 
spectral energy from a broad range of intermediate scales 
into a stationary zonal flow is the mechanism of the zonal 
fiow enhancement during L-H transitions. The result 
shown in Fig.[5]Jb), to large extent confirms this hypoth- 
esis in a radial region close to the transport barrier. Also 
in [3] we presented experimental evidence that the non- 
local spectral transfer of energy from the unstable drift 
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